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I. INTRODUCTION 

It is widely accepted that more than 80% of the matter of our universe is "dark" , a term indicating either that it 
does not interact with photons, or our lack of knowledge of the profound nature of its constituents. Since the seminal 
observations by Oort [1] and Zwicky [2] of the early XX century, there have been several independent evidences that 
some sort of non-baryonic dark matter may exist: the observations of the galactic rotation curves [3] and strong lensing 
data from Hubble space telescope [4] can be accounted for by the single hypothesis that galaxies are surrounded by 
a much larger dark matter halo, which contains most of the mass of the galaxies. This is remarkably consistent with 
the ACDM model of our universe, suggested by the WMAP observations [5], which assumes that the mass energy 
density in our universe is about 30% of the critical energy density, while the observable matter is only 4%, and with 
various observations on clusters of galaxies, which suggest a "mass-to-light" ratio much higher than 1. In addition, 
successful theoretical models and numerical simulations of the formation of structures requires the presence of dark 
matter. 

However, if the existence of dark matter is (almost) commonly accepted, its distribution poses serious challenges. 
In fact, theoretical models and numerical simulations predict, both for clusters of galaxies [6] and for single galaxies 
themselves [7, 8] that the radial density profile is quite steep and peaked at the center of the galaxy. But other 
evidences suggest, on the contrary, that the density profile should be flat, or even shallow [9]. For a comprehensive 
review of the density profile problem (and many others issue related to dark matter), see [10]. 

The possibility that the galactic dynamics of massive test particles may be understood without the need for dark 
matter was also considered in the framework of f(R) modified theories of gravity [11-13]. In particular, the vacuum 
gravitational field equations in f(R) gravity, in the constant velocity region, and the general form of the metric tensor 
is derived in a closed form was analyzed [12]. The resulting modification of the Einstein- Hilbert Lagrangian is of 
the form R 1+n , with the parameter n expressed in terms of the tangential velocity. Therefore, it was concluded 
that in order to explain the motion of test particles around galaxies only requires very mild deviations from classical 
general relativity, and that modified gravity can explain the galactic dynamics without the need of introducing dark 
matter. In the context of modified gravity, the virial theorem in f(R) gravity was generalized by using the collisionless 
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Boltzmann equation [13]. It was found that supplementary geometric terms in the modified Einstein equation provides 
an effective contribution to the gravitational energy. Furthermore, the total virial mass was found to be proportional 
to the effective mass associated with the new geometrical term, which may account for the well-known virial theorem 
mass discrepancy in clusters of galaxies. The model considered in [13] predicts that the geometric mass and its effects 
extend beyond the virial radius of the clusters. Thus, it was shown that the f(R) virial theorem can be an efficient 
tool in observationally testing the viability of this class of generalized gravity models. 

It is interesting to note that recently the possibility that dark matter is a mixture of two non-interacting perfect 
fluids, with different four- velocities and thermodynamic parameters was extensively analyzed [14]. By assuming a non- 
relativistic kinetic model for the dark matter particles, the density profile and the tangential velocity of the dark matter 
mixture were obtained by numerically integrating the gravitational field equations. The cosmological implications of 
the model were also briefly considered, and it was shown that the anisotropic two-fluid model isotropizes in the large 
time limit. In fact, this model was further explored in [15], by assuming that the two dark matter components are 
pressureless, non-comoving fluids. For this particular choice of the equations of state the dark matter distribution 
can be described as a single anisotropic fluid, with vanishing tangential pressure, and non-zero radial pressure. The 
general, radial coordinate dependent, functional relationship between the energy density and the radial pressure was 
also determined, and it was shown to differ from a simple barotropic equation of state. 

Presently, it is considered that the upper bound of the dark matter density in the Solar System is pff M w 2 x 10~ 19 
g/cm 3 = 10 5 (GeV/c 2 ) cm~ 3 [16]. The presence of a local dark matter halo could have some influence on the motion 
of the objects of the solar system [17-19]. The problem of the density and distribution of the dark matter in the 
Solar System is of fundamental importance, not only from a pure theoretical point of view, but also for the design of 
the dark matter particle detectors. Recently an analysis of the kinematics of 412 stars at 1V4 kpc from the Galactic 
midplane performed in [20] has derived a local density of dark matter that is an order of magnitude below standard 
expectations, of the order of 1O _3 M pc -3 , or 0.04 GeV/cm 3 . This result was contested in [21], where it was claimed 
that it arises from the invalid assumption that the mean azimuthal velocity of the stellar tracers is independent of the 
Galactocentric radius at all heights. The assumption of constant mean azimuthal velocity is physically implausible, 
since it requires the circular velocity to drop more steeply than allowed by any plausible mass model, with or without 
dark matter, at large heights above the mid-plane [21]. Using the correct approximation that the circular velocity 
curve is flat in the mid-plane, it was found that the data imply a local dark-matter density of 0.3 ± 0.1 Gev/cm 3 [21], 
consistent with the standard estimates. A new method for the determination of the local disk matter and dark halo 
matter density was proposed in [22]. The method assumes only that the disc is locally in dynamical equilibrium, and 
that the 'tilt' term in the Jeans equations is small up to 1 kpc above the plane. By using this approach a local dark 
matter density of Pdm = 0.025l ; 013 M Q /pc 3 (0.95l ; 49 GeV/cm 3 ) at 90% confidence level, assuming no correction 
for the non-flatness of the local rotation curve, and pnM = 0-022+q;q*3 M Q /pc 3 (0.85j^5o GeV/cm 3 ), if the correction 
is included. The obtained lower bound for the local dark matter density is larger than the standard adopted value, 
and it is inconsistent with the data obtained by extrapolating rotation curves that assume a spherical halo. 

In this paper we will consider a different approach to the problem of the dark matter density in the Solar System. If 
the dark matter density profile follows (roughly) the baryonic matter distribution (for example, because of gravitational 
interaction), it is plausible to assume that there can be specific features at a much lower scale with respect to the 
galactic one. In particular, due to some accretion processes, there can be an excess of dark matter around compact 
objects such as the Sun. Numerical simulation, of course, do not have the resolution to test this hypothesis, but some 
speculations have been made. A brief review of this particular issue is given in [23]. Our aim is to study how classical 
tests of general relativity, i.e., the precession of the perihelion of the planet Mercury, the deflection of light rays and 
the delay of radio signals passing close to the Sun, are modified by the presence of a local, subdominant dark matter 
energy density distribution. The Solar System tests are very powerful tools for testing different gravitational models, 
as well as the theoretical extensions of General Relativity, and they have been applied recently to a variety of contexts 
[24-26]. 

The present paper is organized as follows. In Section II, we derive an approximate expression for the space-time 
metric of the Solar System in the presence of dark matter. In Section III we use this expression, and the formalism 
of Section A to evaluate the corrections due to the presence of dark matter to the relevant parameters in the Solar 
System gravitational tests. Finally in Section IV, we discuss our results and draw our conclusions. In the Appendix 
we describe the general formalism used to evaluate the observational parameters of physical interest. 

II. METRIC OF A STATIC DARK MATTER HALO SURROUNDING THE SUN 



Throughout this work, we assume a general metric of the form 

ds 2 = A(r)c 2 dt 2 - B(r)dr 2 - rW , 



(1) 
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which will be used to obtain the general expressions for the classical tests of general relativity, namely, for the 
precession of the perihelion of Mercury, for the deflection of light rays passing close to the Sun and for the radar echo 
delay observations. For the Schwarzschild metric, giving the spherically symmetric vacuum solution of the Einstein 
gravitational field equations, the functions A(r) = A (r) and B(r) = B (r) arc given by 



(2) 



where rs = 2GM/c 2 is the Schwarzschild radius, with M the mass of the central object. 

We begin this Section by deriving an approximate expression for the metric of a space-time filled by static dark 
matter around a star, by using an ansatz of the form (1). We assume that the dark matter source is described by a 
perfect fluid stress-energy tensor 



T — 



/ \ P( r ) 
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where p and p are related by the equation of state p = jpc 2 , 7 = constant (later on, we will set 7 = according to the 
general equation of state of the matter source) . The Einstein field equations reduce to two independent differential 
equations (we set 87rG/c 4 = k) 



d_ 

dr 



B{r) 

-rA'{r)+A(r) [B(r) - 1] 
r 2 A(r)B(r) 



Kc 2 r 2 p(r), 
—np(r). 



(4) 
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We assume that energy density and pressure of dark matter are very small, p = p\ and p = pi, and hence the 
presence of dark matter can be regarded as a small perturbation of the vacuum state. Therefore we can expand the 
metric around the Schwarzschild solution as 



A(r) = A (r) + A^r), 
B(r) =B (r)+Bi(r), 



Ax(r) < Ao(r), 
Bi(r) « B (r), 



(6) 
(7) 



with A a and B given by Eq. (2). Then, after some algebraic manipulation, using the leading order equation of 
motion and neglecting higher order terms, Eq.(4) becomes 



Bi(r) 



dr [' B 2 (r) 

By defining the effective mass function (see, for example, [27]) 



= kc r pi(r). 



Jo 



m(r) = Air r 2 p(r 2 )dr, 



from Eq. (8) we obtain the function B\(r) as 



B 1 {r)=B 2 {r) 



rs TOi(r) 



(8) 
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where we have expressed the gravitational coupling constant in terms of the mass and of the Schwarzschild radius. 
By using the definition of the effective mass and the equation of state of the matter, we can relate p and p with the 
derivative of m(r). Then, after some algebra and using again leading order equations of motion and solutions, Eq. (5) 
becomes 



J_rs 

M r 



Jmi(r) + y-^ 



TOi(r) 



(11) 



where the prime denotes a derivative with respect to r. 

Now, using the conservation equation to get another independent equation relating the metric and the effective 
mass would be an unnecessary effort, since the barotropic assumption is too constrictive to get a realistic energy 
distribution, and the resulting equations of motion are too difficult to handle. We will instead use, for the dark 
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matter energy density, an effective distribution, inspired by the Navarro-Frenk- White profile for the energy density of 
the galactic halos [19] 

Pl (r) = PBM (—) , (12) 

where pDM and tdm are constant arbitrary parameters related to the mean density and the dimension of the star 
halo, and A is assumed to be a constant. We set 7 = 0, as usual for pressureless matter sources. Then the effective 
mass function and the B\ metric function can be easily evaluated, and are provided by the following relationships 

4^ - - - 3 - X 



m ( r ) = o T^dmPdm (13) 

3- A V dm . 



2-A 



Bl[r) —xT^W^^r- {^) ■ (14) 

A closed expression for A\ can not be found. However, the general solution of Eq. (11) is given by 

M r) = ^^ r ™g^ (l-H) fJ- ^^ dr. (15) 
3 - A r DM M \ r J J r DM (1 - r s /r) 2 

Though it is reasonable to assume A > 0, there is no need in principle to set another limit to A, since our 
approximation holds only outside the Sun, and a matter distribution like (12) can be trusted up to a cutoff radius (for 
example the virial radius in the NFW distribution). Nevertheless, the analysis we are proposing makes sense only if 
the contribution to the metric around the Sun from dark matter is actually localized in the Solar System. Thus, we 
do expect the corrections to the metric functions to vanish as r — > +00, and so we will assume for simplicity A > 2. 

In the next Section, we will use the approximated metric we have deduced to see how the classical tests of general 
relativity constrain the form of the dark matter distribution. 

III. CLASSICAL TESTS OF GENERAL RELATIVITY IN THE DARK MATTER HALO 

In Appendix A, we briefly review the formalism used in studying the classical tests of general relativity, which we 
have included for self-consistency and self-completeness. The formalism is completely general, and can be used to 
treat any static spherically symmetric metric [25, 26]. 

In this Section, we will use the metric given by Eqs. (14) - (15) to study how the results of the classical tests of 
General Relativity (perihelion precession, light deflection and radar echo delay) are modified in the Solar System due 
to the presence of the dark matter. Let us first write down the metric functions in terms of the variable u = 1/r), 

Mu) = o 71 ^77 (l-r S «) / rg-du, (16) 



3 - A r D M M J udm (1 - r s u) 2 

and 

d /. \ ^ 1 r s [dmPdm / u \ A ~ 2 
lW 3-A(l-r s ^r DM M Udm ) ' 1 ' 

respectively. 

A. Perihelion precession 

In order to obtain the change in the perihelion precession due to the presence of dark matter we need to evaluate 
the function G(u) = G (u) + Gi(u), where Go is the standard GR result, given in (A13), and G\ is the first order 
perturbation generated from the modification of the metric functions, which can be written as 

Gi = % (" 2 + h) ^ {AaBi + AiBo) ■ (18) 



■5 



Plugging into this equation the expressions of the metric functions (16)-(17), we arrive at 



Gi = 



47r r s rl M p DM J 2 
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(20) 



respectively. In the small velocity limit, ds ~ cdt, so that E rts c 2 A i rts c (higher order would be subleading in our 
expansion), and we can write 



F = 
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3 — A r DM M 
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(21) 



Using for L and uq the values given in Eqs. (A12)-(A22) respectively, the perihelion precession angle can be written 



as 
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(22) 



The observed value of the perihelion precession for the planet mercury is 8<j>o = 43.11 ± 0.21 arcsec per century 
[28], while the predicted value from GR is 6(j>GR = 42.94 arcsec per century. Assuming that the entire discrepancy 
A5(j) = 0. 17 ± 0.21 is due to the presence of dark matter, we obtain the following constraint for the mean density and 
the radius of the dark matter subhalo 



ruM 
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(23) 



where Tm and Tg are, respectively, the periods of revolution of Mercury and of the Earth. 



B. Light deflection 

To obtain the correction to the deflection of light given by the presence of dark matter around the Sun, we must 
solve the general equation (A20). By defining u = u + u\ + u 2l where w is the straight line solution and u\ is the 
GR result (A27), we can write the equation as 

f^.+u 2 = Q 1 (u ), (24) 

with the functions Q\{u) and P\{u) obtained from the first order corrections to the metric functions due to dark 
matter 1 . For the function Pi we obtain 

P l = is"" 2 - ^Jl ( A ° B 1 + > ( 25 ) 



1 In the right hand side of Eq. (24) we have dropped the term Qo(ui), since it is of order (r 3 /i?) 3 , and can thus be neglected 
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so that, by inserting the solutions (16)-(17), the functions Pi and Qi are given by 
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respectively. 

Keeping only the leading order in the rs expansion, we can again set E 2 /c 2 rts 1 and L w R. Then, inserting the 
expression of u in (27), we can write Eq. (24) as 



d 2 u 2 



U2 



The solution of this equation is 

2W> 3 — A tdm M \ R ) 
The condition u(n/2 + e) = now yields 

__£_ rs_ 2tt r s r DM p DM 1 /r DM \ 
i? i? 2 3 - A r DM M R \ R J 
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(30) 



(31) 



is the Euler's T function [35]. 

Since we have assumed A > 2, we can discard the term e A_1 in the previous equation, so that we are left with a 
simple linear equation. Finally, the deflection angle can be evaluated as 



AGM 
c 2 R 
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n long baseline rad 

ld = S&ld' (1 + ^ld), with Ald < 0.0002 ± 0.0008 arcsec. Thus, assuming, as usual, that all the discrepancy is 

(33) 



The best available data on light deflection by the Sun come from long baseline radio interferometry [29] , which gives 

h (GR) 
LD 

due to the dark matter correction, and taking R = R Q , we obtain the following constraint 

8^ 3/2 HdmPdm r(A/2 — 1/2) /r DM x ' ' 
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C. Radar echo delay 

To obtain the expression of the radar echo delay in the presence of the dark matter sub halo, we have to evaluate 
the integral (A29) with the corrected metric functions. Since we can set 

(34) 
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by neglecting higher order terms we obtain 
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where 2-^1(^,6,0,2;) is the Hypergeometric function [35], defined as the solution of the differential equation 

/- \ d 2 w , , _ N . dw 

z(l - z)— + [c - (a + b + l)z] — -abw = Q 



dz 



(36) 
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Observational constraints on the radar echo delay comes form the frequency shift of radio photons to and from the 
Cassini spacecraft [30]. We have AtRD = At^^(l + Are), with Arr, — (1-1 ± 1-2) x 10 -5 . Therefore, assuming as 
usual that the discrepancy is completely due to the presence of dark matter, we obtain the constraint 
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The next Section will be devoted to some comments on the above results, and we will draw our conclusions. 



IV. COMMENTS AND CONCLUSIONS 



In the previous Section we have evaluated the possible influence that the presence of Sun-bound dark matter could 
have on the classical tests of general relativity. The results we obtained are summarized in Figs. 1 and 2. 

As we can sec, the most constraining test is always the study of the perihelion precession. The results strongly 
depends on the extension of the dark matter halo (Fig. 1). If the distribution is extremely concentrated around the 
Sun, the effect become more and more negligible. On the contrary, for halos extending roughly to the orbit of the 
Earth, the constraint on the mean dark density becomes comparable with the results obtained by other authors [17- 
19]. For halos as big as the planetary solar system, the constraint on the mean dark matter energy becomes tighter. 
On the contrary, at least for what concerns the perihelion precession, the A-dependence is quite mild (Fig. 2). Notice 
that for very large values of tdm the bounds obtained becomes comparable with what is expected from the mean 
energy density distribution of the galaxy [31]. This is a consequence of the choice of our hypothesis for the energy 
density distribution (12), which do not constraint the total amount of dark matter around the Sun. Nevertheless, we 
can argue that, if some sort of overdensity can be localized around the Sun, a distribution packed close to the Sun 
seems to be less in tension with observations. 

Some considerations on how the analysis we have performed can be improved seems appropriate. Of course, 
assuming that the dark matter distribution is static is certainly not true. Probably a model of rotating dark matter 
could improve the constraints, but the metric would be quite complicated. In fact, no exact non-empty spherically 
symmetric rotating solution is known. As a future work, it could be possible to perform a similar analysis, using 
the Kerr geometry as the background expansion. For what concerns the profile energy density, we believe that the 
expression we have used could capture the essential features of a more realistic energy density distribution, especially 
because we are mainly interested in the dimension and the mean density of such a distribution. From this point of 
view, we think that our results are robust. A different but still interesting issue would be to take into account the 
possibility of dark matter "streams" of extragalactic objects recently captured by our galaxy, which could give some 
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FIG. 1: Plots of the constraints on the parameters roM versus pdm of the dark matter distribution p\ of Eq. (12) obtained 
from the perihelion precession (solid line), the deflection of light (dashed line) and the radar echo delay (dotted line). The 
values of A have been set to A = 2.2 for plot (a) and to A = 2.8 for plot (b). The radius is measured in astronomical units, 
while the energy density is measured in GeV/cm 3 . (as a reference for the reader, we remind that the mean value of the galactic 
dark matter halo is of order 1 GeV/cm 3 ) 

strong local overdensities, and which could be characterized by a definite velocity pattern. Maybe it could be possible 
to model such an effect, for example via an Einstein cluster approach [32]. This will also be the object of future 
investigations. 

In conclusion, our analysis suggest that indirect gravitational effect originated by dark matter can be interesting 
even at small scales, and can be a useful cross-check for direct detection experiments. 
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FIG. 2: Plots of the constraints on the parameters pDM versus A of the dark matter distribution pi of Eq. (12) obtained from 
the perihelion precession (solid line), the deflection of light (dashed line) and the radar echo delay (dotted line). The value of 
tdm have been set to tdm — 1AU. The energy density is measured in GeV/cm 3 (as a reference for the reader, we remind that 
the mean value of the galactic dark matter halo is of order 1 GeV/cm 3 ) 
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Appendix A: Classical tests of General Relativity in an arbitrary spherically symmetric static space-time 

In this Appendix, we briefly review the formalism used in studying the classical tests of general relativity. This 
formalism is completely general, and can be used to treat any static spherically symmetric metric [25, 26]. The 
Appendix follows closely the line of the analogous Section of [26], and we include it for self-completeness and self- 
consistency. The general metric of the form 

ds 2 = A{r)c 2 dt 2 - B{r)dr 2 - rW , (Al) 

which will be used to obtain a general expressions for the classical tests of general relativity. For the Schwarzschild 
metric the functions A and B are given by Eqs. (2). 

1. Perihelion precession 

The motion of a test particle in the gravitational field described by the metric given by Eq. (Al) can be derived 
from the variational principle 



S J ^A(r)c 2 i 2 - B(r)r 2 - r 2 (e 2 + sin 2 6»0 2 ) ds = 0, (A2) 

where the dot denotes d/ ds. It is easy to see that the orbit must be planar, and hence we can set 9 = n/2 without any 
loss of generality. Therefore we will use <j) as the angular coordinate. Since neither t nor <f) appear explicitly in Eq. 
(A2), their conjugate momenta yield constants of motion (related respectively with energy and angular momentum 
conservation) 



A(r)c 2 i = E, r 2 4> = L. 



(A3) 
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From (A2), using the constants of motion (A3), changing the variable u = 1/r and expressing the affine derivative 
as an angular derivative d/ds = Lu 2 d/d<p we arrive at the equation 

du\ 2 2 B{u) - 1 2 E 2 1 1 1 , 

T4>) +u =_ ^T U + ^aRbR-^bR = g(u) - (A4) 

By taking the derivative of the previous equation with respect to <j> we find 

dft 

where 



- ns + u = FM, (A5) 



A circular orbit u = u is given by the root of the equation u = F (u )- Any deviation 5 = u — u from a circular 
orbit must satisfy the equation 



AAV 



dF 

du 



5 = O {S 2 ) , (A7) 
which is obtained by substituting u = u + S in Eq. (A5). Therefore, to first order in S, the trajectory is given by 



'-♦■-(f-CfL.H' (A8) 

where 5o and (3 are constants of integration. The angles of the perihelia of the orbit are the angles for which r is 
minimum and hence u or 8 is maximum. Therefore, the variation of the orbital angle from one perihelion to the next 
is 

(A9) 



\J~^ ( du )ii 



1 - cr 



The parameter a defined by the above equation represents the perihelion advance, giving the rate of the perihelion 
change. As the planet advances through <p radians in its orbit, its perihelion advances through <j(f> radians. From Eq. 
(A9), a is given by 



or, for small {dF/du) u=Ual by 



2 V du 



For a complete rotation we have (j> w 2n(l + a) , and the advance of the perihelion is 5<fi — <fi ~ 2tt ~ 2na. In order to 
be able to perform effective calculations of the perihelion precession we need to know the expression of L as a function 
of the orbit parameters. If the planet is moving on a Keplerian ellipse with semi-axis a, and eccentricity of the orbit 
is e, then [25, 26] 

, GMa(l-e 2 ) 
L = y - 2 '-. A12 

cr 

As an example of the application of the present formalism, and for future references, we consider the precession of 
the perihelion of a planet in the Schwarzschild geometry, where we have 

, . 2GM , 1 (E 2 \ 2GM , k , 

G{u) = __ u 3 + (A13) 

, GM 2 GM fk . 

F(u) = 3—u 2 + ^. (A14) 
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The radius of the circular orbit u is obtained as the solution of the quadratic algebraic equation 



GM , GM /A N 

M ° = 3 ~ u o + ^Jj2> ( A15 ) 



with the only physical solution given by 



Therefore 



1±J1- 12G 2 M 2 /c i L 2 GM lk s 

^ = 6GMJ? " *»■ (A16) 



,'dF\ QttGM , A . 



which is the standard general relativistic result [33]. 



2. Deflection of light 

In the absence of external forces a photon follows a null geodesic, ds 2 = 0. The afhne parameter along the photon's 
path can be taken as an arbitrary quantity, and we denote again by a dot the derivatives with respect to the arbitrary 
affinc parameter. There are two constants of motion, the energy E and the angular momentum L, given again by eqs. 
(A3). 

The equation of motion of the photon is 

*+wf?-wf'- (Al8) 

which, with the use of the constants of motion, the change of variable r = 1/u and the use of the conservation 
equations to eliminate the derivative with respect to the afhne parameter leads to 



du\ 2 , B(u) - 1 , E 2 1 . . , KS 

-1 +u 2 = ^J r —u 2 + ^———=P{u). (A19) 



where 



J<pJ B{u) c 2 L 2 A(u)B(u) 

By taking the derivative of the previous equation with respect to <f) we find 

d 2 u 

— +u = Q{u) 1 (A20) 

*> = ^>. ,A 21 , 

In the lowest approximation, in which the term of the right hand side of the equation (A20) is neglected, the solution 
is a straight line, 

cos d> 

«=-^, (A22) 

where R is the distance of the closest approach to the mass. In the next approximation Eq. (A22) is used on the 
right-hand side of Eq. (A20), to give a second order linear inhomogeneous equation of the form 

d 2 u - / cos</>\ 

w +u = Q hr)- (A23) 

with a general solution given by u = u (4>). The light ray comes in from infinity at the asymptotic angle <f> — — (tt/2 + e) 
and goes out to infinity at an asymptotic angle <fi = n/2 + e. The angle e is obtained as a solution of the equation 
u (it/2 + e) — 0, and the total deflection angle of the light ray is 5 = 2e. 
In the case of the Schwarzschild metric we have 

P(u) = (2GM/c 2 )u 3 , (A24) 
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and 

Q(u) = (3GM/c 2 ) u 2 , (A25) 
respectively. In the lowest approximation order from Eqs. (A22) and (A23) we obtain the second order linear equation 



d 2 u 3GM 2 , 3GM .„ „ ,, , A . 

W + U = C ° S * = (1 + C ° S 20) ' (A26) 



with the general solution given by 



cos 6 3GM / 1 „ , , 

+ ' ( A27 ) 



i? 2c 2 i? 2 V 3 



By substituting = 7r/2 + e, m = into Eq. (A27) we obtain the general relativistic result for the light deflection 
angle [33] 

,A 28 , 



3. Radar echo delay 



The time travel of a light signal changes in the presence of a gravitational held. The time difference can be cast as 



l)T = T-T t) = -l iJ^l-l\dy. (A29) 



A(r) 



where we can set r — y/y 2 + R 2 , with R being the distance of closest approach. 
In the case of the Schwarzschild metric we have 



V A(r) c 2 r 
By assuming that the sources are far enough, R 2 jl\ <C 1 and R 2 /l 2 <C 1, and we have [34 



c 3 J- h ^y 2 + R 2 c 3 R 2 
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